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Abstract — In the recent years, a spin-valve was developed with 
a current-confined-path structure consisting of a non-oxide-layer 
(NOL). We analyze magnetic structures of the current-confined- 
path in the nano-oxide layer sandwiched between ferromagnetic 
electrodes and clarify the dependency of the magnetic structure 
on the shape and size of the current-confined-path. Our results 
of stiffness energy density and thermal stability shows that we 
should fabricate a CPP spin-valve with narrow current-confined- 
paths with large aspect ratio and strong constriction in order to 
enhance the MR ratio. 

Index Terms — CPP-GMR, current-confined-path, micromag- 
netic simulation, domain wall 

GEOMETRICALLY confined domain walls of nanometer 
size contacts have attracted enormous attention of re- 
searchers because of its possibility in applications to future 
spintronic devices |ffl- fl6l . In the last decade much effort 
has been devoted to study magnetoresistance (MR) of domain 
walls in atomic wires. One of main obstacles to applications of 
the atomic wires is mechanical stability. For example, Chopra 
et al. 0") reported that a Co atomic wire that shows 300% MR 
ratio is stable for periods of up to as short as 2-3 mins. 

Recently, Fuke et al. fabricated a spin-valve of magnetic 
nanomultilayers with a current-confined-path structure consist- 
ing of a nano-oxide-layer (NOL) with a lot of fine holes filled 
with ferromagnetic metal 04), |fT31 . The holes of an NOL are 
very stable and its typical size is as small as a few nanometers. 
They showed that MR ratio of the spin-valve increases with 
increasing a value of resistance area product, which means 
that the MR ratio due to the geometrically confined domain 
walls is enhanced for the narrow metallic channels. In order 
to understand these experimental results, it is important to 
study magnetic structures of the domain walls confined in the 
NOL. However, it is very difficult to observe the magnetic 
structure of the domain walls directly by using a scanning 
probe microscope or optical techniques because the domain 
walls are confined in the NOL. 

The micromagnetic simulation technique provides a pow- 
erful tool for studying magnetic properties of such geomet- 
rically confined domain walls l6l- |[L3l . Magnetic structures 
of the confined domain walls are investigated in atomic wire 
bridges with crystalline anisotropy by several groups l6l- lfl0l . 
However the effect of the crystalline anisotropy on magnetic 
structures should be small for the NOL spin-valve fl4l . Ifl5l 
where the size of the domain walls is of the order of a few 
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nanometers. Although a few simulations for contacts without 
crystalline anisotropy have been reported in |[TT1 - |[T3"1 . a certain 
careful simulation where the real shape of the current-confine- 
path in the NOL spin-valve lfT4l . Ifl5l is taken into account is 
needed to design the highly sensitive magnetoresistive devices. 

In the present paper, we study shape dependence of mag- 
netic structures of current-confined-paths in the NOL by using 
the micromagnetic simulation. We evaluate stiffness energy 
density of the magnetic structures in the NOL sandwiched 
between two magnetic layers of which magnetizations are 
aligned antiparallel to each other. We can estimate the ef- 
fective thickness of the geomehically confined domain wall 
from the calculated stiffness energy density. We show that a 
narrow current-confined-path with large aspect ratio and strong 
constriction is required in order to enhance the MR ratio. 

We employ a spin model on the simple cubic lattice, which 
is defined by the Hamiltonian, 



H 



-K d J2$i ■ J drD{xi - r) ■ S(r) 

i 

-iv-„^(l-(e l -S ? 4 ) 2 )+^^-^ 



(1) 



The local values at the i-th site, B%, Xi, e, and hi denote an unit 
vector representing a classical Heisenberg spin, a coordinate, 
an unit vector aligned in a direction of crystalline anisotropy 
and an applied magnetic field, respectively. The first term in 
the right hand side of (Q]i is the exchange stiffness energy 
between spins at nearest neighbor sites. The exchange coupling 
constant is denoted by J. The value of J is related to exchange 
stiffness constant, A, and a lattice constant, a, as J = 2aA. 

The second term is the dipole-dipole interaction energy. The 
tensor function, T>{x), is expressed as, 

1 1 — 2>es ® eg 



V(x) 



An 



where eg = x/\x\. The dipole-dipole coupling constant, Kd, is 
equivalent to a 3 /xoM^/2, where M s and (j,q are the saturation 
magnetization and vacuum permeability, respectively. In order 
to deal with the dipole-dipole interaction in the NOL spin- 
valve, we adopt a finite element - boundary element (FEM- 
BEM) hybrid method |"P71 . In the evaluation of the stray 
field coming from the dipole-dipole interaction, spins at lattice 
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Fig. 1. The finite element mesh of the current-confined-path in the NOL 
spin-valve is shown. The shape of the current-confined-path is deformed to 
a hyperboloid parametrized by two ratios, d/h and d c /d, where d, d c and 
h, denote the diameter of the center, that of the bottom and height, of the 
cun'ent-confined-path, respectively. 



Fig. 2. The magnetization configurations of the Neel-wall like ground state 
for d/h = 0.8 is shown in the panel (a). The magnetization configurations of 
the Bloch-wall like ground state for d/h = 0.8 is shown in the panel (a). The 
diagram of magnetic structures of the ground states is shown in the panel (c). 
The Neel-wall like , Bloch-wall like structures are indicated by the diamonds, 
cubes and triangles. 



points are extrapolated to that in the whole region in the NOL, 
S(f). 

The third term is the crystalline anisotropic energy. We ne- 
glect this term in the present calculation because the character- 
istic length determined by the crystalline anisotropy is a few- 
tens of nanometers and therefore the effect of the crystalline 
anisotropy on the magnetic structure of the geometrically 
confined domain wall of a few nanometer size must be very 
small. 

The forth term is the Zeeman energy due to the applied 
magnetic field. We also neglect the Zeeman energy in the 
present calculation and the direction of the magnetizations in 
the top and bottom electrodes are determined by the boundary 
conditions. The boundary conditions we use are the following: 
Sf = —1 on top and side surfaces of the top electrode and 
Sf = 1 on bottom and side surfaces of the bottom electrode. 
On the other surface the boundary condition is free. 

The system we consider is schematically shown in Fig. 03 
The system is divided to a 3-dimensional lattice consisting of 
hexahedral finite elements. The number of elements in the top 
electrode is taken to be 39 x 39 x 3 = 4563. The bottom 
electrode is divided in the same way as in the top electrode. 
The number of elements in the z-direction in the contact is 
taken to be 9. The total number of the elements, which depends 
on the ratio d/h, is of the order of 10 4 . The shape of the 
contact is assumed to be a hyperboloid parametrized by two 
ratios, where d and d c are the diameters, respectively, at the 
center and the bottom of the contact, and h is the height 
of the contact. The characteristic length determined by the 
competition between the exchange interaction H ex and the 
dipole-dipole interaction Hd is the exchange length defined 

by, 

Z ex = \/4-, (2) 
V K-d 

which is of the order of a few nanometers for conventional 
ferromagnetic metals 11811 . We assume that Z cx =2.8nm cor- 
responding to the materials used in 1141 . The unit length of 
the system, a, is taken to be y/Kd/Jl ex and the height of the 
NOL is h = 2nm (~ OJliex). 

The classical Heisenberg spin system is embedded in the 
whole region as shown in Fig. Q] The spins obeys the Landau- 
Lifshitz-Gilbert equation, 

d ?t 7 7t r dH * dH , 
dt l + a 2 l dS, dSi 

where 7 and a are the gyromagnetic ratio and the Gilbert 
damping constant, respectively. The calculated ground states 
by the present simulation do not depend on their values. We 
simulates relaxation of the spins from various initial states 
by numerically solving the equation at zero temperature to 
obtain the ground state. The time integration is continued 



Fig. 3. Concentration of stiffness energy at d c /d= 0.5 d/h= 0.8 (a), d c /d= 
0.5 d/h= 3.9 (b), d c /d= 1.0 d/h= 0.8 (c) and d c /d= 1.0 d/h= 3.9 (d). 

until a space- and time-averaged value of the torque exerted 
to the spins is of the order of 10~ 8 7J/(1 + a 2 ). In the 
present simulation, further time-integration does not produce 
additional change of magnetic structures within the reachable 
time scale by our computers. The obtained lowest energy state 
by simulations started from three initial states is regarded as 
the ground state at each given set of parameters. 

In the present simulation, by changing two ratios, d/h and 
d c /d, we investigate the systems in the range of 0.8 < d/h < 
4.0 and for d c /d = 0.25, 0.5, 0.75 and 1.0. The obtained ground 
states are the Neel- and Bloch-wall like states as shown in 
Figs. |2 a) and|2|b), respectively. 

The magnetic structure diagram is shown in Fig. |2£c). At 
small values of d/h, the ground states are the Neel-wall 
like state. On the other hand the ground states observed at 
large values of d/h are the Bloch-wall like state. The similar 
results were obtained for atomic wires with a crystalline 
anisotropy by Coey et al. J6|. The result can be understood 
as follows: Our results show that the magnetic structures 
depend almost only on z-component of coordinates. Thus the 
exchange energy does not depend on whether the Bloch-wall 
or Neel-wall like structures J3) and we can discuss the energy 
of the states only by the dipole-dipole interaction energy. The 
dipole-dipole interaction energy on the surface of the contact, 
which decreases with increasing the aspect ratio, d/h, gives a 
dominant contribution to the total energy of the Bloch-wall like 
state. On the other hand, the dipole-dipole interaction energy 
in the interior of the contact, which increases with increasing 
the aspect ratio, d/h, gives a dominant contribution to the 
total energy of the Neel-wall like state. Therefore the Bloch- 
wall (Neel-wall) like state is preferred for the wide (narrow) 
contact as show in Fig. 2 (c). As we increase the ratio d c /d, the 
dipole-dipole interaction on the surface of the contact increases 
but that in the interior of the contact decreases. Therefore the 
value of the aspect ratio d/h which determines the boundary 
between the Bloch- and Neel-like wall states is a decreasing 
function of the ratio d c /d. 

Figs. Oa)-[3jd) show the dependence of the stiffness energy 
density on the size and shape of the contact. The stiffness 
energy density is defined as 

E^±ldr(vS { r)f, (4) 

~ A P dz\V0{z)\\ (5) 
where 6(z) denotes azimuthal angle for the Bloch-wall like 
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state and polar angle for the Neel-wall like state. Vi is total 
volume of elements near the i-th site. The last expression in (O 
is justified for much small values of Kd/ J. In our simulation, 
value of Kd/ J about 0.006 and thereby the approximation is 
also justified. 

The stiffness energy densities of the narrow contacts are 
shown in Figs. (3ja) and 0b). The stiffness energy densities 
are almost confined in the NOL and its spacial distribution 
can be controlled by changing the shape of the contact. 
Comparing the results shown in Figs. 0a) and Ob), one can 
see that the stiffness energy density of d c /d = 0.5 is much 
more concentrated at the center of the system than that of 
d c /d = 1.0. Since MR ratio increases with decreasing the 
thickness of the domain wall fl6l . we can expect a larger value 
of the MR ratio of the contact with d c /d = 0.5 than that with 
d c /d = 1.0. 

The stiffness energy densities of the wide contacts are 
shown in Figs. [3jc) and [3jd). The stiffness energy densities 
stick out from the NOL and therefore the effective thickness 
of the domain wall is thicker than that for the system with the 
narrow contacts. As shown in Fig. 0c), a strong concentration 
of the stiffness energy density appears around the surface of 
the contact with d c /d=0.5. Although a strong concentration 
of the stiffness energy density means that we have a thin 
domain wall around the surface, we cannot expect the large 
MR ratio because the other region with small stiffness energy 
density gives a small contribution to the MR ratio and acts as 
a parasitic resistance. 

Finally, we discuss the stability of the ground states. The 
bistability of the Bloch and Neel-like structures is observed 
for almost all of the region of the parameters except for d c /d 
=0.25 and d/h < 3. Thus thermal instability is expected 
in those as discussed in J6), 0. In our results, the energy 
difference between the Bloch- and Neel-wall like states is 
of the order of a few Kelvins at most. Although in wide 
contacts the energy difference between the Bloch- and Neel- 
wall like states increases with increasing system volume. If 
h is set at the typical experimental value of 2 nm fl4l . 
d should be a few tens nanometers at least in order that 
the energy difference reaches room temperatures. However, 
as discussed before, we can not expect large MR ratios for 
such wide contacts. In the Neel-wall like states, we do not 
observe the excited Bloch state for d c /d =0.25 and d/h < 3.3 
within our calculation accuracy. Therefore the narrow contacts 
with large aspect ratio and strong constriction has stable 
magnetization. The narrow contacts are a potential candidate 
of highly sensitive magnetoresistive devices because of their 
high magnetoresistance ratio as mentioned before. 

In summary, we analyze magnetic structures of a current- 
confined-path in a nano-oxide layer sandwiched between 
ferromagnetic electrodes and clarify the dependency of the 
magnetic structure on the shape and size of the current- 
confined-path. Our results of stiffness energy density and ther- 
mal stability shows that we should fabricate narrow current- 
confined-path with large aspect ratio and strong constriction 
in order to enhance the MR ratio. 
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